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Abstract. 

The quasar B 2005+403 located behind the Cygnus region in our Galaxy, is a suitable object for studying the interplay be- 
tween propagation effects, which are extrinsic to the source and source intrinsic variability. On the basis of VLBI experiments 
performed at 1.6, 5, 8, 15, 22, and 43 GHz during the years 1992 - 2003 and parallel multi-frequency monitoring of the total 
flux density (between 5-37 GHz), we investigated the variability of total flux density and source structure. Below 8 GHz, the 
point-like VLBI source is affected by scatter-broadening of the turbulent interstellar medium (ISM), which is located along 
the line of sight and likely associated with the Cygnus region. We present and discuss the measured frequency dependence 
of the source size, which shows a power-law with slope of -1.91 ± 0.05. From the measured scattering angle at 1 GHz of 
9 = 77. 1 ± 4.0 mas a scattering measure of SM= 0.43 ± 0.04 irT 20 ' 3 kpc is derived, consistent with the general properties of 
the ISM in this direction. The decreasing effect of angular broadening towards higher frequencies allows to study the inter- 
nal structure of the source. Above 8 GHz new VLBI observations reveal a one-sided slightly south-bending core-jet structure, 
with stationary and apparent superluminally moving jet components. The observed velocities range from 6.3 - 16.8 c. The jet 
components move on non-ballistic trajectories. In AGN, total flux density variations are often related to the emergence of new 
VLBI components. However, during almost eleven years no new component was ejected in B 2005+403. A striking feature in 
the flux density variability is a trough observed at 5 - 37 GHz and between 1996 and 2001. This trough is more pronounced at 
higher frequencies, where it shows a prominent flux density decrease of ~ 50 - 60 % on a time scale of 2 - 3 yrs. The trough 
can be explained as a blending effect of decreasing and increasing jet component fluxes. Dense in time sampled flux density 
monitoring observations with the 100 m Effelsberg telescope reval intra-day variability (IDV) at 1.6 GHz with a modulation 
index of m = 1.0 %. At 5 GHz less pronounced variations are seen (m = 0.5 %). This and a relatively short variability time scale 
of ~ 0.1 days imply a second, less dense or turbulent scattering screen at nearby (few to hundred parsec) distance. 
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1. Introduction 

The propagation of radio waves through the ionized interstel- 
lar medium causes several effects, such as Faraday rotation 
and depolarization of polarized emission, dispersion of pul- 
sar signals, scatter broadening of compact radio sources and 
intensity fluctuations caused by diffractive (DISS) and refrac- 
tive (RISS) interstellar scintillation (ISS). The effect of scatter 
broadening plays an important role in the understanding of the 
compact object located in the center of our galaxy Sgr A* (e.g . 
iDoeleman et al.l200lllKrichbaum et alll998bHLo et al.ll 998). 
It is however also important for other galactic regions on the 
sky, e.g. for the Cygnus region and the quasar B 2005+403 lo- 
cated behind. The detailed study of scatter broadening in com- 
pact radio sources could lead to a better understanding of the 
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interstellar medium. It also bears the potential to reveal fine 
details of the intrinsic source structure, particularly at longer 
wavelengths, where direct imaging often cannot provide suffi- 
ciently high angular resolution. 

The other prominent propagation effect addr essed in this 
paper is the so called Intra-D ay Variability (IDV, I Witzel etafl 
1986; iHeeschen et all 1 19871) . and the question if and how 
it is related to the interstellar scintillation. Recent IDV- 
surveys show that a large fraction (up to 30%) of all com- 
pact flat spectrum radio sources show this effect, which 
is characterized by variability amplitudes of up to 20 - 
30 % and variability times cales ranging front less than one 
hour t o several days (e.g. iKraus et alj 120031 IHeesc hen et alJ 
Il987l iKedziora-Chudczer et all l200ll lLovell et alJ l2003l)~ lf 
interpreted via source intrinsic incoherent emission pro- 
cesses, such short variability timescales imply - via the light 
travel time argument - apparent brightness temperatures of 
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Tn = 10 16 - 19 K, far in excess of the inverse-Compton limit 
jKellermann & Paulinv-Tothlll969l) . With the assumption of 
relativistic Doppler-boosting the brightness temperatures can 
be reduced, which however, would imply uncomfortably large 
Doppler boosting factors of 6 - 20 - 200. 

Alternately, the IDV in the radio-bands is also explained 
extrinsically via scintillation of radio waves in th e turbu- 
lent interstellar medium (ISM) of our Galaxy (e.g. iRickett 
1990; Dennet-Thorpe & de Bruvn 2000). The main problem 
in this interpretation is that it lacks the explanation of ob- 
served ra dio-optical broad-band c orrelation s in at least some 
source s fouirrenba ch et alJ dl99 lh . see also IWaener & Witzell 
l l 1995b V The recent detection of diffractive ISS in the source 
J 1819+3845 leads to micro-arcsecond source sizes and bright- 
ness temperatures of ~ 10 14 K, which again requires Dopple r 
boosting factors of 6 - 100 JMacquart & de Bruvr]l2005l) . 
It is therefore likely that the IDV phenomenon invokes both, 
a co mbinatio n of source intrinsic and extrinsic effects (e.g. 
iKrichbaum et al 120021: loian et al.l2002l) . 

In this paper we present new measurements for the scatter 
broadened quasar B 2005+403, combining Very Long Baseline 
Interferometry (VLBI) and flux density variability measure- 
ments obtained at various frequencies over the last decade. 
B 2005+403 is a flat spectru m quasar (5s oh?. = 2.5 - 3.5 Jy, 
tto.3/5GHz = 0. 3 with S ~ y^lBeckeretal.lll99ll) at a redshift 
of z — 1.736 (B oksenberg et alJll976l) . It is located close to 
the Galactic plane at / = 76.82°, b = A.29° behind the Cygnus 
super-bubble region. Earlier studies showed that interstellar 
scattering affects the VLBI image of the source, causin g an- 
gular broadening at frequencies below 5 GHz llFev et all 1989: 
iMutel & Lestradell990tlDesai & Fevfcool . The high value of 
the scattering measure (SM = f Cl(s)ds '), de rived for the 
line of sight of B 2005+403 bv lFev et al.1 Jl989i> - reflects the 
strong influence of interstellar medium. 

In order to further investigate the scatter broadening and 
the amount of scintillation in B 2005+403, we combined all 
available flux density monitoring data with our VLBI obser- 
vations. The wide frequency coverage of our data facilitates a 
study of the interplay of the frequency dependent scattering ef- 
fects, which dominates below 8 GHz and the source intrinsic 
variability, which is best seen at and above 15 GHz. At these 
higher frequencies B 2005+403 has not been intensively stud- 
ied with VLBI before. The VLBI data cover a time range of 1 1 
yrs (from 1992 - 2003) and allow us to study and monitor the 
structural evolution of VLBI-core and pc-scale jet. The scatter 
broadening observed in the VLBI images at the lower frequen- 
cies is measured and gives further constraints to the properties 
of the intervening ISM. Additional parameters for the ISM are 
obtained from the IDV monitoring of the source, performed 
with the Effelsberg 100 m telescope. 

The paper is organized as follows. In Sect. 2, we present 
details of the VLBI observations and corresponding data re- 
duction and introduce the long-term trends of the flux density 
variability. In Sect. 3 we discuss the source structure, the kine- 
matics of the inner jet components and the possible location of 

1 The Scattering Measure is the path integral over the coefficient C' N 
of the electron density fluctuation wavenumber spectrum. 



the nucleus (VLBI core). In Sect. 3.4, we focus on the results 
of the low frequency VLBI observations and discuss the scatter 
broadening. In section 3.5 we discuss the long-term variability 
of B 2005+403 and relate it to changes of the VLBI structure. 
In Sect. 3.6 we discuss the implications from the observed in- 
traday variability. Sect. 4 gives a concluding summary. 

The following cosmological param- 

eters are used throughout this paper: 
ff = 71 kms- 1 Mpc" 1 , Q m = 0.27 and Q vac = 0.73. 

2. Data and data reduction 

2. 1. VLBI data 

We analyzed 17 VLBI datasets of B 2005+403 from 14 dif- 
ferent observing epochs obtained during 1992 - 2003. These 
observations were performed with different VLBI networks 
(EVN, VLBA, Global) at various frequencies ranging from 
1 .6 GHz to 43 GHz. In our analysis we also included several 
data sets from the 2 cm-survey (see iKellermann et alJ 12004) 
and one observation o f the source in the MOJAVE survey 
jLister & HomaiJl2005l) . In Table □ we summarize the VLBI 
experiments. 

After correlation at the VLBI correlators in Socorro 
(NRAO) or Bonn (MPIfR), the VLBI data were calibrated 
and fringe-fitted using the standard procedures within AIPS 
(Astronomical Image Processing System) software package. 
The post-processing of the data invoked the usual steps of edit- 
ing, phase and amplitude self-calibration and imaging. These 
tasks were performed within AIPS and the Caltech DIFMAP 
packages. For previously published datasets (from the 2 cm- 
survey and the MOJAVE survey and the datasets from epoch 
2003.04), we reanalyzed the data. Starting from pre-calibrated 
visibilities, we redid the self-calibration, imaging and Gaussian 
model fitting. For comparison of the results from this reanalysis 
w ith the original maps in t he 2 cm Survey we refer the reader 
to IKellermann et alJ J2004 . 

In three VLBI experiments the source was also observed in 
full polarization. Column 4 of Table^summarizes the recorded 
polarization. In the polarization analysis, the feed leakage 
terms for the antennas (D-terms) were determined in the usual 
way using the AIPS task "LPCAL" (Lepp anen et alJll995l) . At 
15 GHz, we could use the University of Michigan (UMRAO) 
database to calibrate the orientation of the electric vector po- 
sition angle (EVPA). We verified that total and correlated flux 
density (in total intensity and in polarization) matched within 
10% accuracy, excluding a large miscalibration of the EVPA 
owing to overseen polarized large scale structure. At 22 and 
43 GHz, we had no possibility to calibrate the EVPA. At these 
frequencies, the absolute orientation of the EVPA is therefore 
unknown. 

In Fig.^ we present the results from the VLBI imaging, 
showing the CLEAN maps obtained at 15 GHz, 22 GHz, and 
43 GHz. Figure^also shows the polarization images with the 
electric vectors superimposed. 

To facilitate measurement of structural changes, the VLBI 
maps were parameterized in the usual manner using Gaussian 
components. We fitted Gaussian components to the calibrated 
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visibility data using the DIFMAP program. The resulting mod- 
elfits are shown in a separate table which is available online. 
The table contains the following information: Col. 1 gives the 
observing epoch, frequency and VLBI array, Col. 2 identifies 
the VLBI components by letters and numbers, Col. 3 gives the 
flux density of each VLBI component (S), Col. 4 its radial dis- 
tance r relative to component CI, Col. 5 the position angle 
(6), Col. 6 the FWHM major axis (a), Col. 7 the ratio of mi- 
nor to major axis (b/a), Col. 8 the position angle of the major 
axis (0), and Col. 9 the reduced x 1 to characterize the good- 
ness of the modelfit. The positional information (Columns 4 
and 5) are relative to the westernmost component, denoted as 
CI. It is further assumed that the position of CI is stationary 
with time. The measurements errors listed in the table were es- 
timated from the internal scatter obtained during the process of 
imag ing and model-fitting a nd is described in more detail e.g. 
in Krichb aum et alJ Jl99 8a). During the process of model fit- 
ting we reduced the number of independent variables by fitting 
only circular Gaussian components. For those cases where cir- 
cular Gaussian model fits gave not satisfactory results and pro- 
duced unacceptable large reduced^ 2 values, elliptical Gaussian 
components were fitted. Most of the 22 GHz, and 43 GHz data 
were fitted best with elliptical Gaussians. 

In Figure 13 we show the VLBI images resulting from the 
Gaussian model fitting. The circles seen in this figure denote for 
the positions and sizes of the individual model fit components. 
The contour lines result from the convolution of the modelfit 
components with the observing beam and the addition of the 
underlying (noise limited) residuals. 

2.2. Total flux density monitoring 

B 2005+403 is included in the flux density monitoring 
programs performed at the University of Michigan Radio 
Astronomy Observatory (UMRAO) in the USA and at the 
Metsahovi Radio Observatory in Finland. This allows to study 
the long-term flux density variability. In Figure|7]we show the 
lightcurves at 5 GHz, 8 GHz and 15 GHz jAller et all 19991) and 
at 22 GHz and 37 GHz (Terasrant a et all 19 98: Teras ranta et alJ 
2004). The covered time range is 1990 - 2004. A more de- 
tailed discussion of the variability and its possible relation to 
the structural changes seen in the VLBI maps is given in Sect. 
3.5. 

In the context of ongoing systematic studies of Intraday 
Variable radio sources, B 2005+403 was observed with the 
Effelsberg 100 m radio telescope of the Max-Planck-Institut fur 
Radioastronomie. During two observing runs especially dedi- 
cated to search for IDV in sources located in the Cygnus re- 
gion (the observed sources are those studied bv lFev et al.ll989l 
iDesai & Fevl l200ll) . we measured the flux density variabil- 
ity of B 2005+403 at 6 cm (4.85 GHz) and 18 cm (1.67 GHz) 
with the Effelsberg 100 m radio telescope. The observations 
were performed on November 30 and December 6, 2003. 
The flux densities were measured adopting the method of re- 
peated and subsequently averaged cross-scans. The measure- 
ment method and the dat a analysis is described in more detail 
e.g. in lKraus et al.l II2003I) . The duty cycle of the measurement 



resulted in 1 -2 flux density measurements per source and hour. 
For calibration purposes, several steep spectrum sources were 
included and observed with the same duty cycle as the target 
source B 2005+403. From these additional sources, the source 
B 2021+614 showed the least variability and therefore could be 
used as secondary calibrator to monitor and correct for residual 
short time-scale fluctuations of the antenna gain. The absolute 
flux density scale was determined from repeated observations 
of the primary calibrato rs NGC7027, 3C 48 and 3C286 and 
the flux density scale of dBaars et al.ll977HOtt et alll994l) . 

3. Results 

3. 1. Kinematics 

In Fig.[T]we show the VLBI CLEAN maps obtained at 15 GHz, 
22 GHz, and 43 GHz. The corresponding maps of the Gaussian 
modelfits are presented in Fig. [2] 

The maps of B 2005+403 show an east-west oriented core- 
jet structure with several embedded components. The central 
1 mas region is best described by three Gaussian components 
CI, C2, and C3. The relative alignment of these 3 components 
seems to vary over the time spanned by the observations. In 
most images, the relative alignment of CI - C3 suggests a slight 
bending to the north. Beyond 1 mas core separation and ori- 
ented more to the south-east, a faint and partially resolved re- 
gion of diffuse and extended emission is visible. It can be fitted 
by one or two Gaussian components of larger extend (C4 and 
C5). 

In Fig. [5] we plot the relative separation from CI for the 
components C2 to C5 versus time. Owing to its inverted spec- 
trum, compactness and variability, we assume that CI is the 
VLBI core and that its position is suitable for use as a station- 
ary reference point (Sect. 3.2 for a discussion). In the figure, 
we also superimpose a straight line resulting from a linear fit to 
the motion of each component. The slope of each line measures 
the angular separation rate. These rates and the corresponding 
apparent velocities are summarized in Table|2] 

With the exception of component C2, which appears to be 
stationary within the measurement errors, the components C3 
to C5 move with apparent superluminal velocities. For these 
three components we observe a systematic increase to their ap- 
parent velocities from 6.3 c for component C3, to 9.9 c for C4 
and 16.8c for C5. In the 2 cm survey jKellermann et alJ20 04). 
B 2005+403 is characterized by a two component model, with 
the secondary component moving at an apparent speed of 
12.3 + 3.0c. Based on its core-separation, this component may 
be identified with a blending of components C4 and C5 in our 
analysis. 

The minimum Lorentz factor y m ; n of the jet could be 
estimated from the maximum observed component speed. 
Adopting the usual equations for 




and the inclination angle of the jet, which leads to maximum 
apparent speed 

lp 7mm = COt^jSapp, 



4 



K. E. Gabanyi et al.: High frequency VLBI observations of the scatter broadened quasar B 2005+403 



Table 1. The epochs, frequencies and arrays of the VLBI observations described in this paper. In Col. 4 we indicated the polar- 
ization mode of the observation. In Col. 5 references for previously published data are given. 



Epoch 


Frequency 


Instrument 


Pol. 


Reference 


1992.44 


22 GHz 


EVN (8 stations) 


LL 


this paper 


1994.17 


22 GHz 


VLBA (8 stations) + VLA + EVN (4 stations) 


LL 


this paper 


1995.27 


15 GHz 


VLBA (2 cm Survey) 


RR 


Kellermann et al. (2004) 


1995.96 


15 GHz 


VLBA (2 cm Survey) 


LL 


Kellermann et al. (2004) 


1996.73 


43 GHz 


VLBA+EB 


full 


this paper 


1996.73 


22 GHz 


VLBA+EB 


full 


this paper 


1996.73 


15 GHz 


VLBA+EB 


full 


this paper 


1996.82 


5 GHz 


FVN (8 stationsi 


LL 


this nanpr 

1111.3 I'Cll. 1 ^! 


1996.83 


8 GHz 


FVN C8 stationsi 


RR 


this nanpr 

1111.3 I'Cll. 1 ^! 


1997.19 


15 GHz 


VLBA (2 cm Survey) 


LL 


this paper 


1998.14 


1.6 GHz 


EVN (9 stations) 


LL 


this paper 


1999.01 


15 GHz 


VLBA (2 cm Survey) 


LL 


Kellermann et al. (2004) 


2001.17 


15 GHz 


VLBA (2 cm Survey) 


LL 


Kellermann et al. (2004) 


2001.98 


15 GHz 


VLBA (2 cm Survey) 


LL 


Kovalev et al. (2005) 


2003.04 


22 GHz 


VLBA 


full 


Bach (2004) 


2003.04 


15 GHz 


VLBA 


full 


Bach (2004) 


2003.16 


15 GHz 


VLBA (MOJAVE) 


full 


Lister & Homan (2005); Kovalev et al. (2005) 



one obtains for component C5 y m ; n > 16.8 + 2.3 and i/r rmin = 

3.3 ° + 0.5°. 

Laht eenmaki & Valtaoial \ 19991) calculated the Doppler 
boosting factor for a sample of active galactic nuclei using the 
total flux density monitoring data of the Metsahovi observatory. 
For B 2005+403, they report a variability Doppler boosting fac- 
tor of (5 var = 8.63, which translates to £ var = 10.9, for the cosmo- 
logical parameters used in this paper. We note that because of 
the lig ht travel time argument used bv lLahteenmaki & Valt aoia 
Jl999l) . the derived Doppler-factor places a lower limit to the 
true Doppler-factor. With the proper motions obtained from 
VLBI, we may calculate the viewing angle and the Lorentz fac- 
tor of the jet: 

2/? app 

^ var = arctan — — — — 

Papp ' ^var A 

_ lapp + <C + 1 
7var — ~ c 

zo var 

(Avar < 4.8° + 0.9° and y var > 18.4 + 3.5. 

The apparent acceleration observed for components C2 to 
C5 may either result from an intrinsic acceleration of the jet, or 
from a systematic bending of the jet, or from a combination of 
both effects. As t// 7min and i[> v . dr both are very small, a jet bending 
of less than a few degrees could easily explain the observed in- 
crease of the apparent velocities. To investigate this further, we 
plot in the lower panel of Fig.[3]the apparent variations of the 
position angle of the radius vector for each component. The 
plot indicates that component C2 and C3 move along similar 
and spatially bent trajectories with position angle changes in 
the range of ~ 20 - 40°. Between 1992 and 2003 the position 
angle difference between C2 and C3 increases, indicating still 
similar but spatially offset paths. For components C4 and C5 
the situation is less clear, as much smaller variations of the po- 
sition angles are observed and the scatter in the data is larger. 
We therefore conclude that at least the components C2 and C3 
move on non-ballistical and spatially bent trajectories. 



Table 2. The apparent speed and the corresponding apparent 
proper motion of components C2, C3, C4 and C5. 



Component Apparent speed Apparent proper motion 





mas yr ' 


in units of c 


C2 


(-0.3 ± 4.0) 10" 3 


< 0.31 


C3 


0.082 + 0.006 


6.25 ± 0.46 


C4 


0.14 ±0.03 


9.91 ±2.12 


C5 


0.22 + 0.03 


16.77 ±2.29 



The variation of the flux density of the individual jet com- 
ponents is plotted versus time in Fig. [4] Since most data are 
available at 15 GHz, we connect the flux density measurements 
at this frequency with lines. The variations of the flux density 
of components CI, C2 and C3 appear completely uncorrelated, 
with the flux density of CI mainly decreasing, that of C2 in- 
creasing and the flux density of C3 showing some oscillatory 
behavior with at least two local maxima. In Section 3.5 we will 
discuss the combined flux density variability of the VLBI com- 
ponents with the variations seen in the total flux density of the 
source. 

3.2. The position of the VLBI core 

From our simultaneous multi-frequency VLBI observations the 
identification of the VLBI core and the measurement of in- 
dividual component spectra is facilitated. For those compo- 
nents, where quasi-simultaneous flux density measurements 
were available at different frequencies, we calculated the spec- 
tral indices and show in Table the results. In 1996.73 and 
2003.04, component CI shows between 15 GHz and 22 GHz 
an inverted spectrum. Between 22 and 43 GHz, the spectrum 
steepen. Components C2 and C3, both show an optically-thin 

synchrotron spectrum, with similar spectral indices of 0.7 

consistently measured in both observations. The spectral shape 
of CI can be explained via synchrotron-selfabsorption. This, 



K. E. Gabanyi et al.: High frequency VLBI observations of the scatter broadened quasar B 2005+403 



5 
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1 (b) 
1995.96 
15 GHz 






Kg) 
" 2001.98 
15 GHz 





Hi) 

2003.16 
15 GHz 

1.560 Jy/bcam 
2.34 mjy/bcam 




a 



I® 

-1992.44 
22 GHz 



0.584 Jy/bcam j 
5.84 mjy'/beam 




l(k) 
1994.17 
22 GHz 



0.448 Jy/bcam 
4.48 mjy/bcam 




I 





Fig.l. VLBI CLEAN maps of B 2005+403 at 15 GHz, 22 GHz, and 43 GHz shown at the different observing epochs. Epoch, 
frequency, peak flux density and the lowest (positive) contour are given in the upper left corner of each image. Contours are 
in percent of the peak flux and increase by factors of two. The beamsize is shown in the bottom left corner of each image. 1 
mas length of the superimposed polarization vectors corresponds to 12.5, 25, 10, 25, 31.3, 31.3 mJy/beam in images (h) to (n), 
respectively. (Images were created in AIPS.) 



its flux density variability and the high compactness seen in 
the modelfits indicate, that CI has to be identified with the 
synchrotron-selfabsorbed jet base, which commonly is called 
the VLBI core. 

We note that despite its optically thin spectrum, compo- 
nent C2 shows the most pronounced flux density variability 



(see Figure|4}. This is somewhat atypical, as pronounced flux 
density variability is usually seen in jet components with much 
flatter spectra. The apparent flux density increase of AS isghz = 
1 .3 Jy in Af var = 2.16yrs leads via the light travel time argu- 
ment to an apparent brightness temperature of = 2.2-10 12 K, 
which apparently violates the inverse Compton limit. The cal- 
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10.86 mjy/beam 



Fig. 2. Maps of Gaussian modelfits for B 2005+403. Circles characterize positions and sizes of the individual modelfit com- 
ponents. Contours result from the convolution with the observing beam. Epoch, frequency, peak flux density and the lowest 
(positive) contour level are given in the upper left corner of each image. Contours correspond to the modelfits, they are in percent 
of the peak flux and increase by factors of two. As in Fig.^ the VLBI observations are at 15 GHz, 22 GHz, and 43 GHz. (Images 
were created in DIFMAP.) 
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ml 1 1 ' — — i 1 1 1 1 1 1 ' 1 

1992 1994 1996 1998 2000 2002 2004 



Epoch [year] 

Fig. 3. Top panel: core separation as a function of time for 
components: C2 (circles), C3 (triangles), C4 (diamonds), C5 
(left triangles). Filling of a symbol denotes the frequency of 
the observation: open for 22 GHz, filled for 15 GHz, larger 
and striped for 43 GHz. Solid lines indicate linear regres- 
sion to reliable identifications (C2 and C3). The dashed-dotted 
lines show regression to tentatively suggested identifications. 
Bottom panel: components position angle as a function of time. 
Same symbols and lines are used as above. Lines are used just 
to guide the eye. 



culation lead to a lower limit for the Doppler-factor of 5 > 3.6, 
which would be required to reduce the observed brightness 
temperature to 10 12 K. It is therefore very likely that the vari- 
ability is mainly caused by differential Doppler boosting in 
combination with motion along a curved path. The relatively 
large variability and the apparent stationarity of C2 strongly 
indicate motion more or less directly towards the observer, i.e. 
an angle to the line of sight close to zero. A similar interpre- 
tation to explain apparent station arity has been ma de also for 
other sources, i.e. for 4C 39.25 jAlberdi et all2000l) . 
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Fig. 4. Flux density evolution of individual jet components CI, 
C2 and C3 is shown. CI is assumed to be the core. Open 
grey symbols denote observations at 22 GHz, filled black sym- 
bols at 15 GHz and open diamonds represent the observation at 
43 GHz. 

Table 3. Spectral indices of the components CI, C2, and C3 

(S ~ v +a ) 



Epoch 




CI 


C2 




C3 


1996.73 


< 

< 


= 1.34 
= -0.35 


af s = -0.70 


af 5 


= -0.70 


2003.04 




= 1.46 


af 5 = -0.65 




= -0.67 



On the basis of their spectral shapes with a steep spectral 
index of -0.7, we rule out that components C2 or C3 can be 
identified with the VLBI core. An implication of this is, that 
B 2005+403 is a one-sided core-jet source and that it exhibits 
no counter-jet. The lack of a counter-jet is usually explained as 
a consequence of Doppler-boosting, consistent with our obser- 
vations. 

The spectral shape of CI indicates the existence of a spec- 
tral maximum near 22 GHz (Table 0. If this is identified 
with the synchrotron-turnover of a homogeneous synchrotron 
self-absorbed component, we can estimate the magnetic field 
strength using the measured size of component CI: 0.34 mas. 

Following (e.g. Marscher 1983), the magnetic field strength 

is: 

B = 10 J <5-^-7 -T 

S 2 (l+z) 

where 6 is the Doppler factor, 8 is the component size in units 
of mas, v is the turnover frequency in GHz, S (in Jy) is the flux 
density at the turnover frequency and b(a) is a spectral inde x 
dependent tabulated parameter (see Table 1 in Mars cherll98l . 
Here we adopt a = -0.35. This yields for C\B=\.15- 10~ 4 T. 
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Table 4. Polarization characteristics of B 2005+403 obtained 
from the VLBI maps. Col. 1 gives the epoch, Col. 2 the fre- 
quency of the observation. Col 3 gives the name of the compo- 
nent. Col. 4 gives the polarized flux density in mJy, Col. 5 the 
polarization angle in degrees. The EVPA is not calibrated the 
22 GHz and 43 GHz datasets; the given values are arbitrary. 



Epoch 


v [GHz] 


ID 


P [mly] 


xO 


2003.04 


15 


C2 


66+14 


131 ±3 


2003.04 


22 


C2 


74+14 


131 ±3 


2003.16 


15 


C2 


67 ± 14 


121 ±6 


1996.73 


15 


C3 


42 ±5 


75 ±8 


1996.73 


22 


C3 


35 ±4 


49 ±8 


1996.73 


43 


C3 


13 ±3 


90 ±4 


2003.04 


15 


C3 


38 ±4 


90 ±4 


2003.04 


22 


C3 


21 ±3 


90 ±4 


2003.16 


15 


C3 


28 ±5 


90 ±5 




3.3. Polarization 

In Fig. [0 the polarization vectors are superimposed on the 
VLBI maps for those epochs where we have polarization data. 
The lengths of the vectors represent the polarized flux density. 
The polarized flux and the polarization angles derived from 
the VLBI observations are displayed in Table 0] Owing to the 
lack of absolute polarization angle calibration at 22 GHz, and 
43 GHz, the position angles of the electric vectors are arbitrary 
at these frequencies. The measured polarization characteristics 
at 15 GHz are in agreement with the UMRAO single dish mea- 
surements. 

Component C2 was unpolarized in 1996.73, however this 
was the most prominent polarized feature in 2003.04 and 
2003. 16, with comparable polarized flux density at 15 GHz and 
22 GHz. The polarized flux density of C3 decreased through all 
the epochs and it also decreased by increasing frequency. 

The polarization angles of C2 did not change signifi- 
cantly with time. The polarization angle of C3 changed by » 
15° and 40°, at 15 GHz, and 22 GHz respectively, from 1996.73 
to 2003.04. In 2003.04 and 2003.16 it remained at the same 
value, roughly parallel with the jet direction implying that the 
magnetic field is perpendicular to the jet, as one expects from 
an optically thin component. 

The absolute angle of polarization is unknown at 22 GHz 
and 43 GHz however the angles in all epochs for both compo- 
nents are not significantly different from the calibrated EVPA 
at 15 GHz (one exception is 1996.73 at 22 GHz). This is consis- 
te nt with the low Rotation Measure 2 (RM) values determined 
bv lZavala & TavloJ ( l2003l) . 

3.4. Low frequency VLBI observations 

The central part (denoted as Ml in Table [TJ of B 2005+403 
cannot be resolved at 1.6 GHz, 5 GHz, and 8 GHz. However, 
at 5 GHz and 8 GHz two and one extended features appear in 



100 50 -50 

Right Ascension (mas) 



100 



2 RM is an integration of magnetic field strength multiplied by the 
electron density along the line of sight from the source to the observer: 
RM = 812 rc e B||ii/radnr 2 , where B\\ is the parallel component of 
the magnetic field in mG. 



Fig. 5. Clean map obtained in EVN observation (at 1 .64 GHz) 
in 1998.14. The peak flux density is 0.464 Jy/beam. The con- 
tour levels are 4.64 mJyfteam x(-3.2, 3.2, 6.4, 12.8, 25.6, 51.2). 
The beam FWHM is 13 mas, in contrast to the apparent angular 
size of the source of 36 mas. 



the jet direction, towards south-east (M2 and M3). M2 can be 
tentatively identified with C5 in the 15 GHz and 22 GHz data. 

Our low frequency VLBI observations at 1.6 GHz, 5 GHz 
and 8 GHz confirm the effect of ang ul ar broaden i ng rep orted 
earlier bv IMutel & Lestradd dl990l) . iFev et alJ |i989) and 
besai & Fevl feOOll) . In Figure we show 1.6 GHz map ob- 
tained with the EVN in 1998.14. The image illustrates the ef- 
fect of scatter broadening: the measured source size (36 mas) is 
three times larger than the size of the observing beam (13 mas). 
The possible elongation seen at this map in the position angle 
of ~ 37° may indicate anisotropy in the scattering. 

Similar to lFev et al.l Jl989l) . we analyze the measured total 
angular source size as a function of the observing frequency. 
For this, we fitted the visibility data of B 2005+403 at all avail- 
able frequencies with one Gaussian component. In Figure [6] 
we plot the FWHM size of this component versus frequency, 
using our data obtained at 1.6 GHz, 5 GHz, 8 GHz, 15 GHz, 
22 GHz, and 43 GHz and including the previous measurements 
from the literature dFev et al.lll989t IMutel & Lestradd 1 1990t 
Desa i & FevlEool . the corresponding values are tabulated in 
Table |5] The figure clearly shows a deviation from a simple 
power-law. In order to quantify the scatter broadening, which 
dominates at lower frequencies, we fitted a power law to the 
data in the frequency range of 0.67 - 8 GHz. In this range we 
obtain the following result: 

0= (77.1 + 4.0) • (y/l GHz)- (1 - 91±005) 

The dashed line in Figure |6] shows this fit. We note that 
the slope does not change significantly, if the data points at 
8 GHz are excluded from the fit. In this case one obtains: 
6 = (11.1 ± 4.0) • (v/1 GHz)" 1 - 92 * 06 , In both cases the slope 



K. E. Gabanyi et al.: High frequency VLBI observations of the scatter broadened quasar B 2005+403 



9 



of the power-law differs significantly from the slope of -2.2 
one would expect for Kolmogorov type density fluctuations in 
the ISM. 

Above 8 GHz the source is significantly larger than the ex- 
trapolated scattering size. This indicates that towards higher 
frequencies scattering effects are less dominant and that the 
intrinsic structure of the sources shines through. The differ- 
ences between the extrapolated scattering size and the mea- 
sured source size are as follows: A#i5ghz = (0.4 ± 0.1) mas, 
A6»22GHz = (0.4 ± 0.1) mas, and A6> 43 ghz = (0.4 + 0.1)mas. 
For the three size measurements above 8 GHz, we fit an- 
other power law, which to first order can be used to charac- 
terize the frequency dependence of the intrinsic source size: 
Bint = (1.7 +0.4) ■ (v/1 GHz)- (a31±0 07) mas. The solid line in Fig. 
[6] shows this fit. With this curve we are able to obtain an upper 
limit of the intrinsic source size also at lower frequencies, were 
direct size measurements with VLBI are not possible. 

Because of its heavily scattered line of sight, B 2005+403 
was incorporated by Cordes and Lazio in their model of 
describin g the d istribution of the Galactic free electrons 
JCordes & Laziol Eo02). To account for the regions of intense 
scattering, they included "clumps", which are regions of en- 
hanced electron density and/or el ectron density fl uctuation. 
Adopting the scattering measure of lFev etafl 1^989) thickness 
d » 18 pc located at a distance ofL* 2.35 kpc was derived. In 
their model, a Kolmogorov spectrum of the interstellar turbu- 
lence was assumed. 

At a given frequency, the size of the scattering disk can be 

estimated via the deconvolution formula: 6* scat = - # 2 nt . 

For the scattering si ze at 1 GHz, we obtain Q\ gh? = (77.1 ± 
4.0) mas. Following Tavl or & Cordel ll 19931) and assuming 
Kolmogorov turbulence, we obtain for the scattering measure: 



SM 



/ "scat \ 1 

1 128 mas/ ' V<i 



1/3 
GHz' 



, 128 mas; 

The derived value of SM = (0.43 + 0.04)irT 20/3 kpc is in 
goo d agreement and consistent with the previous measurement 
of lFev etalJ fi989), however now has a smaller uncertainty. 
Our slightly lower scattering measure implies either a some- 
what reduced electron d ensity fluctuation than that given in the 
Cordes & Lazio (2002) model or a shorter path length through 
the screen. For the el ectron density fluctuation as given in 
ICordes & Laziol (|2*002), we formally obtain a path length of 
1 1 . 7pc, instead of the 18 pc me ntioned above. 

Spangler & Re ynolds! 111990) reports on the measurements 
of Ha intensity towards the Cygnus region. They studied lines 
of sight towards extragalactic radio sources which are affected 
by interstellar scattering. In the direction of B 2005+403, the 
reported Ha intensity is 130 Rayleigh, which corresponds to 
an emission measure of EM = 293.8 cirT 6 pc (assuming an 
electron temperature of 8000 K). The EM corresponds to the 
path integral of the squared electron density along the line of 
sight. For a thin screen it can be approximated by: EM = 
n 2 ■ d. Assuming that the probed HII region cau ses the scatter 
broadening using equation (16) and (18) from lCordes & Laziol 
(2002) we calculated from the EM and SM values the outer 
scale to be Iq = 25.5 AU, and the fractional variance of the 



1 10 100 

Frequency [GHz] 

Fig. 6. The measured angular size plotted versus observing fre- 
quency. The dashed line represents a power law fit to the data 
in the range of 0.6 GHz to 8 GHz. The slope of the line is 
— 1.91 +0.05. The solid line represents a power law with a slope 
of -0.31, fitted to the data points at the high frequencies in the 
range of 15 GHz to 43 GHz. Circles denote data from the liter- 
ature, triangles data from this paper. 



electron density to be 6.9. Using the thickness of 11.7pc the 
electron density can be calculated: n e = 5.0cm -3 . 

The exponent of the power-law of the angular broadening 
(-1.91 ±0.05) indicates a wavenumber spectrum of the electron 
density turbulence, which is flatter than the usually assumed 
Kolmogorov type turbulence with spectral slope of -2.2. This 
can be produced in a medium consisting of random superpo- 
sitions of disc o ntinuities such as a c ollection of shock fronts 
dRicketti ri990). Ro mani et ail |l986) tabulate numerical esti- 
mates of scintillation parameters for different power-law spec- 
tra of electron density fluctuations. W e used the equation w ith a 
power law index of -2 from Table 1 of Ro mani et alJ (fl986) that 
describes most adequately our angular broadening observation. 
From this we obtain C\d = 15.16rrr 20 ^ 3 kpc. With a thickness 
of the screen in the range of d — 12 - 18 pc we obtain a range 
for C 2 N = (842... 1266) rrT 20 / 3 . Ho wever, we note that this is not 
an entirely consistent approach, as Cordes & Lazio (2002) used 
a Kolmogorov spectrum. 

From the power-law fitted to the source size at high fre- 
quencies, an upper limit for the intrinsic source size and hence 
lower limit for the brightness temperature can be calculated. 
At 1 .6 GHz, 5 GHz, and 8 GHz the intrinsic source sizes are 
f3 int < (1.5 ± 0.6) mas, int < (1.0 ± 0.3) mas, and 6 mt < 
(0.9 ± 0.2) mas, respectively. The corresponding lower limits 
to the brightness temperatures are 7b > (0.6 ± 0.5) ■ 10 12 K, 
T B > (0.14 ± 0.09) • 10 12 K, and T B > (0.08 ± 0.04) • 10 12 K. 
These numbers are in accordanc e with typical b rightness tem- 
peratures measured in VLBI (e.g. lKovalev et al . 2005) and nei- 
ther strongly violate the inverse Compton limit nor indicate ex- 
cessive Doppler-boosting. 
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Table 5. The measured angular sizes of B 2005+403 at differ- 
ent frequencies. Col. 1 shows the observing frequency, Col. 2 
the observing epoch, Col. 3 the major axis size of the fitted 
Gaussian component. In Col. 4 references for previously pub- 
lished data are given. 



v [GHz] 


Fnnrh 




T3 p-ppi-pripp 


0.61 


Oct 1986 


225 ± 58 


Fev et al. (1989) 


1 fx 

1.0 


1 QQS 1 A 

1 yyo. 14 


Zy.o ± V.J 


this paper 


1.66 


Mar 1986 


26.4 ±4.7 


Fevetal. (19891 


1.7 


Jan 1997 


31 ±0.8 


Desai & Fev (2001) 


2.3 


Jan 1997 


17.6 ±0.5 


Desai & Fev (20011 


4.99 


Oct 1985 


3.5 ±0.4 


Fev et al. (19891 


5 


1996.82 


3.5 ±0.1 


this paper 


5 


Jan 1997 


3.4 ±0.1 


Desai & Fev (20011 


8 


1996.83 


1.6 ± 0.1 


this paper 


15 


1996.73 


0.8 ±0.1 


this paper 


22 


1996.73 


0.6 ±0.1 


this paper 


43 


1996.73 


0.5 ± 0.04 


this paper 



3.5. Long-term flux density variations 

In Figure we show the radio light-curves of B 2005+403 at 
5 GHz, 8 GHz, 15 GHz, 22 GHz, and 37 GHz during the time 
interval covered by our VLBI observations. The data between 
5 GHz and 15 GHz are from the Michigan (UMRAO) moni- 
toring program, the data at 22 GHz and 37 GHz are from the 
Metsahovi monitoring. 

The data clearly show the longterm variability of the 
source. A feature of particular interest is the apparent decrease 
of the flux density observed between 1992 and 1999. Between 
1999 and 2001, the flux increased again. This flux density 
'trough' is best seen and more pronounced at the higher fre- 
quencies. At 5 GHz, the decline begins at 1996.5 and the flux 
density drops by 16 % reaching its lowest value around 1999.4. 
At 8 GHz, a flux density decrease of 30 % is observed during 
approximately the same time interval. At 15 GHz, the flux den- 
sity drops by 47 % from 1996.6 till the beginning of 1999. At 
the two highest observed frequencies (22 GHz and 37 GHz), the 
decrease is even more pronounced: 54 % and 60 % respectively. 
This systematic frequency dependence of the variability ampli- 
tude is accompanied by a systematic and frequency dependent 
shift of the time of the flux density minimum. It is obvious that 
the flux density minimum and the subsequent rise of the flux 
density appear earlier at the higher frequencies. In Fig. we il- 
lustrate this frequency dependence by a solid line, which is not 
a fit and only meant to guide the eye. 

To quantify the variations at the different frequencies we 
use the modulation index, defined as follows: 

m = 100— 

(S) 

%, were cr denotes for the rms variation and {S ) for the mean 
flux density. The mean flux density and the corresponding vari- 
ability indices are summarized in Table [6] It is seen that at 
22 GHz the modulation index is largest. 

In a number of AGN, short time flux density troughs are ex- 
plained via an occultation by ISM clouds moving through the 
line of sight. This are the so called "extreme scattering events" 



Table 6. The mean flux density, the standard deviation and 
the modulation indices for the flux density variability of 
B 2005+430 at different frequencies between 1990 and 2004. 



Frequency (GHz) 


<5>(Jy) 


cr 


m(%) 


5 


2.86 ± 0.02 


0.217 


7.61 


8 


3.03 ± 0.02 


0.369 


12.15 


15 


2.63 ± 0.02 


0.374 


14.24 


22 


1.95 ±0.03 


0.500 


25.67 


37 


1.84 ±0.03 


0.363 


19.81 



(ESE, see lFiedler et al.ll994l) . Depending on size, distance and 
relative velocity of the scatterer, timescale be tween a few days, 
months and perhaps even years are possible JCimo et all2002t 
IPohl et alJfl995l) . Since B 2005+403 is located in a region of 
the sky with prominent scattering effects, the interpretation of 
the observed flux density trough via such an extreme scatter- 
ing event provides an attractive possibility for an interpreta- 
tion. However, the frequency dependence and variability time 
scales seen in B 2005+403 are different than those expected 
for an ESE. Typical ESE events show more pronounced and 
more rapid variations towards longer wavelengths, opposite to 
what is observed in B 2005+403. Also the frequency depen- 
dence of the flux density minimum observed in the trough of 
B 2005+403 is not consistent with an ESE event, which owing 
to the geometrical occultation should not show any frequency 
shifts. 

The observing dates of the low frequency VLBI observa- 
tions performed in 1996.8 at 5 GHz and 1998.1 at 1.6 GHz co- 
incide with the times of the decrease of the total flux density. 
The observed angular source sizes at these times, however, are 
not significantly differe nt than the sizes observed previously 
(1985 & 1986) bv lFev etalJ Jl989l) and long before the flux 
density trough. This indicates that the amount of the observed 
angular broadening is persistent and therefore is not related to 
the observed flux density variability. 

In an alternative interpretation we argue that the observed 
flux density trough can be explained by the sum of the flux 
densities of the variable VLBI components. In Figure [7] we 
display the sum of the flux density of the modelfit compo- 
nents C1+C2+C3 from the VLBI observations at 15 GHz and 
22 GHz. Particularly at 15 GHz, were the time coverage of the 
VLBI experiments covers the full time range of the trough, 
there is an excellent agreement between the single dish flux 
density measurements and summed VLBI component fluxes. 
Examining the flux densities of the individual components (see 
Fig-H}, it appears as if the combination of the declining trend 
of CI and the increasing trend of C2 are mainly responsible 
for the shape of the flux density trough. We therefore conclude 
that the peculiar shape of the total flux density lightcurve of 
B 2005+403 is a result of a blending of the independently vary- 
ing flux of the inner jet components. 

3.6. Short time-scale flux density variations 

In November and December 2003 the flux density variability 
of B 2005+403 was monitored with the Effelsberg 100 meter 
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Fig. 7. Light-curves of B 2005+403 at 5 GHz, 8 GHz, 15 GHz, 22 GHz, and 37 GHz from the monitoring data of the University 
of Michigan Radio Astronomy Observatory, and the Metsahovi 14 meter radio telescope. The straight line connects the minima. 
The dark circles and lines denote the sum of flux density of the modelfit components from the VLBI observations at 15 GHz and 
22 GHz. 



radio telescope with high time resolution. In Fig.[8]and Fig.|9| 
we plot the variability curves at 18 cm and 6 cm, respectively. 

In Table we summarize the results from these measure- 
ments. To characterize the variability amplitudes and thei r sig- 
nifican ce we follow the meth ods described in iHeeschen et alJ 
dl987h and lKraus et~ai1 ( l2003l) . Column 2 of the table gives the 
average flux density, col. 3 the rms standard deviation, col. 4 
the modulation index, col. 5 the variability amplitude (defined 

as Y = 3 Jm 2 - m^, see iRraus et al-lEoO^ and col. 6 the re- 
duced x 1 - A value of Y is only printed, if the ^ 2 -test gives a 
higher than 99.9 probability for significant variations, mo repre- 
sents the maximum modulation index of the calibrators. From 
table [7J and Figures [3] and |5] it is obvious that the variability 
amplitude of B 2005+403 decreases with frequency. Formally 



we measure a modulation index of m = 1.01 % at 1.67 GHz 
and of m = 0.45 % at 5 GHz. At this frequency the detection 
of IDV is marginal. We note that the model of refractive inter- 
stellar scintillation in the weak regime predicts a dec rease of 
the m odulation index with increasing frequency (e.g. iRickettl 
Il990l) . This is consistent with our observations. 

The characteristic variability timescales derived from the 
light curves in Figures [8] and grange between f var = 0.09 - 
0.13 days at 18 cm and f var = 0.05 -0.1 2 days at 6 cm. Adopting 
a source intrinsic interpretation for the variability we may apply 

the light travel time argument and derive - via the source size 

I 1 

- an a pparent brightness temperature. Following Wagner et al. 
( 1996) we obtain for the apparent "variability brightness tem- 
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12974.1 12974.2 
J.D.-2440000 



Fig. 8. Light-curve of B 2005+403 at 1.6 GHz in November 
2003 observed with the Effelsberg 100 meter telescope. Filled 
circles denote B 2005+403, open squares (connected with 
dashed line) the calibrators, NGC7027 and B 2021+614. The 
modulation indices of the secondary calibrators are a measure 
of the calibration uncertainty, which is too < 0.25 %. The mod- 
ulation index of B 2005+403 is m = 1.01 %. 



1.02 - 
1.015 
„ 1.01 - 




12981.1 12981.2 
J.D.-2440000 



Fig. 9. Light-curve of B 2005+403 at 5 GHz in December 2003 
observed with the Effelsberg 100 meter telescope. Filled circles 
denote B 2005+403, open squares (connected with dashed line) 
a calibrator, NGC 7027. The modulation index of the calibrator 
is mo = 0.20 %, the modulation index of B 2005+403 is to = 
0.45 %. 

Table 7. The variability parameters of B 2005+403 and the sec- 
ondary calibrators at 1.67 GHz (top) and 4.85 GHz (bottom). 
Col. 2 shows the mean flux density in Jy, col. 3 its standard de- 
viation, col. 4 the modulation index, col. 5 the noise-bias cor- 
rected variability amplitude and col. 6 the reduced^- 2 . 



Source name 


<S>(Jy) 


£7" 


m(%) 


Y(%) 


red.^ 2 


A = 18 cm, m = 0.25% 


B 2005+403 
NGC 7027 
B 2021+614 


2.430 
1.906 
2.193 


0.024 
0.004 
0.006 


1.01 
0.19 
0.25 


2.93 




5.593 
0.197 
0.365 


A = 6 cm, m = 0.20% 


B 2005+403 
NGC 7027 


2.905 
5.489 


0.013 
0.011 


0.45 
0.20 


1.18 



3.107 
0.621 



T B = 4.5 • 10 10 AS 



AD, 



fvar(l +ZY 



K 



where AS is the amplitude of the flux density variations [Jy], 
A [cm] the wavelength, Dl [Mpc] the luminosity distance, and 
? V ar [day] is the variability time-scale. The resulting apparent 
brightness temperatures is in the range of Tb = 10 18 ■ 21 K. This 
can be reduced to the inverse-Compton limit with Doppler- 
factors in the range of a few hundred of up to one thousand. 
Since these values appear unreasonably high and are not con- 
sistent with the results fr om VLB I nor with theoretical expecta- 
tions (e.g. lBegelman et al.ll994h . we probably can rule out that 
source intrinsic effects are responsible for the observed rapid 
variations. 

The source size and the apparent brightness temperature 
can be also determined from the interstellar scintillation model, 
assuming that the main reason for the observed variability is 
the motion of the Earth thro ugh the scintillation pattern (see 
for example iGoodmanll 19971 and references therein). In this 
scenario the variability timescale (in [days]) is given by: f var = 
347.10 scat LtT I with the scattering size 8 sc . dt in [mas], the screen 
distance L in [kpc] and the relative velocity between screen and 
earth v in [km/ si . It is clear that a scree n located at kpc distance, 
as proposed by Cordes & Lazio (2002) is not able to explain the 
observed short variability timescale. With a scattering size (at 
1.67 GHz) of # scat = 30 mas, a screen distance of 2.4 kpc and 
typical galactic velocities u < 220km/s, one would expect to 
see variations on timescales f var > 100 days. 

In order to reproduce a variability time scale of 0. 1 day 
at 1.67 GHz we obtain with the above equation the follow- 
ing constraint for scattering size and screen distance: # scat L < 
6.34 • 10~ 2 . For the measured scattering size of ~ 30 mas, this 
leads to an unreasonable nearby screen distances of < 2pc. 
The only way out of this dilemma is a smaller scattering size. 
Adopting a minimum screen distance of at least 10-20pc as 
required for the int erpretatio n of the ultra fast scintillators (e.g. 
Dennett-Thorpe & de Bruvn 2002, 2003), we would obtain a 
scattering size of scat < 6.3 mas. A source with such a size 
would have an apparent brightness temperature at 1 .6 GHz of 
7b > 10 10 K, fully consistent with our VLBI measurements 
(see Section 3.4). 

An upper limit to the distance of the screen is obtained 
from the restriction of the source size via the inverse Compton 
limit of Jg C = 10 12 K. The requirement that the brightness 
temperature has to be lower than this, leads to a source size 

of 9 = ^1.22 • ^SyV^T^d-Hl + z) > 0.5 mas, where we 
adopted Si.eGHz = 2.4 Jy from Table and for the Doppler 
factor 6 - 10. This leads with the above equations to an 
upper limit of the screen distance of D < 0.58upc or with 
(10 < v < 220)km/s (6 < L < 127)pc. 

We therefore conclude that the observed IDV can be ex- 
plained with a screen, which is located at a distance of L < 
130pc and which can be characterized by a scattering size of 
0.5 < # scat < 6.3 mas. The corresponding scattering measure at 
1.6 GHz then is of order 6.4- 10~ 4 < SM < 4.3 • 10~ 2 , measured 
in the conventional units of SM. These values are considerably 



K. E. Gabanyi et al.: High frequency VLBI observations of the scatter broadened quasar B 2005+403 



13 



lower than the scattering measure of the more distant screen, 
which is responsible for the scatter broadening in B 2005+403. 
We therefore conclude that the observed IDV is not caused by 
this distant screen. Most likely one has to deal with the situation 
of multiple scattering by at least two spatially and physically 
very different plasma screens. In this scenario the first screen 
leads to a significant scatter broadening of the source image, 
which in the second screen is only very weakly scattered due to 
large quen ching effects (fl Mim >> # SC at). This quenched scat- 
tering (e.g. lRickettl dl986l) ') by the second screen could explain 
the relatively low variability amplitudes of < 1 % observed in 
our IDV experiments at 6 and 18 cm. Quantitatively this can be 
verified using equation (20) of Goodman (1997), which relates 
the variablity index, the scattering measure, the effective source 
size and the screen distance. With the parameters from above a 
modulation index of m < 1 % are obtained, in good agreement 
with our observations. 

4. Discussion and summary 

B 2005+403 is located in the Cygnus-region and previously 
was known as one of the most scattered broadened extragalac- 
tic sources in this region of the sky. In this paper we combined 
multi-frequency flux density monitoring and VLBI imaging 
observations. We also used high time resolution flux density 
monitoring observations, in order to search for possible Intra- 
Day Variability. The combination of these complementary data 
yields new and more accurate than previously published esti- 
mates of (i) the source structure and its kinematics and of (ii) 
the properties of the interstellar medium, which is responsible 
for the observed propagation effects. 

High frequency VLBI imaging observations (at 15 - 43 
GHz) spanning 11 years (1992 - 2003) reveal a one sided and 
bent core jet structure, with at least 5 embedded VLBI com- 
ponents. The inner jet components separate from the station- 
ary assumed core CI with apparent superluminal speeds of 6 
- 17 c. The component nearest to the core (C2), however, ap- 
peared stationary (yS app < 0.3). It is remarkable that the flux 
density variation of this stationary and steep spectrum compo- 
nent apparently violates the inverse Compton limit, indicating 
strong Doppler-boosting. This could be explained with a com- 
ponent path, which is oriented at a smaller angle to the line 
of sight than those of the other jet components. A systematic 
increase of the component speeds with increasing core separa- 
tion is observed. The paths of the inner-most jet components is 
curved and their motion is not ballistical, suggesting a spatially 
curved (helical) jet. Beyond 2 mas core-separation, however, 
the components seem to move on linear (ballistical) trajecto- 
ries. Continued future VLBI monitoring will be necessary to 
determine their paths more accurately. 

The VLBI images at longer wavelengths (1.6 - 8 GHz) 
show a scatter broadened nearly point-like source. Combining 
our new size measurements with published data from earlier 
observations allowed to study the frequency dependence of the 
source size in greater detail. Below 8 GHz this dependence is 
best described by a non-Kolmogorov power law with slope of 
-1.9. For the size of the scattering disk at 1 GHz we obtain 
(77.1 + 4.0) mas. Above 8 GHz the measured sizes are increas- 



ingly larger than the prediction from the scattering law, and 
the internal source structure becomes visible. Based on the ob- 
served scatter broadening, several parameters of the scattering 
medium (scattering size, scattering measure, electron density) 
were determined and found to be in accordance with previ- 
ously published estimates (NE2001 model, Cordes & Lazio, 
2002). The NE2001 model places the scattering screen at a 
distance of 2.35 kpc. We however note that towards the di- 
rection of B 2005+403 there is a degree-size extended super- 
nova remnant G78.2+2.1. It exhibits a patchy, inhomogeneous 
struct ure at optical, X-ray and radio wavelengths (B vkov et alJ 
2004 1 and references th erein). Its distance is (1.5 +0.5) kpc 
(Landecker etail [?980 ), still compatible with the scattering 
model. 

We discuss the observed long term flux density variability 
curves (5 - 37 GHz), which during 1996 - 2001 showed a re- 
markable flux density 'trough'. A tentative interpretation via 
a possible 'scattering event', similar to the 'extreme scattering 
events' observed in other sources, could be abandoned on the 
basis of the observed frequency dependence of the variabil- 
ity, which appeared more pronounced at higher frequencies. 
Instead, we show that the 'trough' is well explained by the 
summed flux density variability of the (independently) vary- 
ing VLBI components (CI - C3). In contrast to other sources, 
which often show correlations between jet comp onent ejection 
and flux density variability ('e.g. lSavolainen et alJ2002l and ref- 
erences therein), such a relation is not seen B 2005+403. In this 
source at least some (if not most) of the observed flux density 
variability results from the blending of the evolving jet compo- 
nents. 

Dense in time sampled variability measurements performed 
with the 100 m Effelsberg telescope showed only weak intra- 
day variability in B 2005+403. This was in contrast to naive 
expectations, in which more pronounced variability was an- 
ticipated. The low variability indices (m = 1% at 18 cm, 
m = 0.5 % at 6 cm) and the short variability time scale of 
-0.1 days cannot be explained by scattering effects from the 
kpc-screen, which is responsible for the scatter broadening. 
Instead, another and much closer scatterer is required. Using 
the thin screen approximation for refractive interstellar scintil- 
lation we determined the likely distance of the scattering screen 
to be (6 < L < 128) pc. At 18 cm, its scattering size would lie in 
the range of 0.5 < 6» scat < 6.3 mas. For B 2005+403 one there- 
fore has to deal with the situation that a scattered broadened 
image (kpc-screen) is scattered a second time (pc-screen) and 
that the combination of both effects limits the accuracy of low 
frequency flux density variability measurements to about 0.5- 
1 % and the angular the resolution of VLBI observations to a 
few milliarcseconds. We note that weak refractive interstellar 
scintillation from a very nearby and ionized ISM (which most 
likely covers larger portions of the sky), may therefore also 
limit the calibration accuracy and sensitivity of planned large 
radio telescopes operating and low frequencies, like LOFAR or 
the square-kilometer array (SKA). 
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Table .1. The parameters of the Gaussian modelfits. Column 1 gives the observing epoch, the VLBI array and the frequency, column 2 the 
component identification. Column 3 shows the flux density of each component. In Col. 4, the relative separation to CI and in Col. 5 the position 
angle are listed. Columns 6, 7 and 8 list the major axis, the axial ratio of minor to major axis, and the major axis position angle. The last column 
shows the reduced^ 2 of the model. 
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Epoch, array, frequency 
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